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Abstract 

In this paper, we claim the availability of deterministic noises for stabi¬ 
lization of the origins of dynamical systems, provided that the noises have 
unbounded variations. To achieve the result, we first consider the system 
representations based on rough path analysis; then, we provide the notion 
of asymptotic stability in roughness to analyze the stability for the sys¬ 
tems. In the procedure, we also confirm that the system representations 
include stochastic differential equations; we also found that asymptotic 
stability in roughness is the same property as uniform almost sure asymp¬ 
totic stability provided by Bardi and Cesaroni. After the discussion, we 
confirm that there is a case that deterministic noises are capable of mak¬ 
ing the origin become asymptotically stable in roughness while stochastic 
noises do not achieve the same stabilization results. 


1 Introduction 

Asymptotic stability (AS) is an important property for analyzing the motion 
of dynamical systems. In control theory, the property ensures that trajectories 
of systems converges on the desired steady-states liiiniiii]. In a particular 
case, the convergence is achieved by just adding white noises; the strategy that 
directly aims at this phenomena is said to be stabilization by noise [DiiiiTiiin]- 
The strategy of stabilization by noise provides the origins with asymptotically 
stable in probability; however, the stability is weaker than AS because it does 
not ensure the existence of invariant sets, which are proper closed subsets of 
Euclidean spaces. The stochastic stability property “almost the same” as AS is 
uniform almost sure asymptotic stability (UASAS) [3111]; all the sublevel sets of 
the related Lyapunov functions are invariant sets with probability one. However, 
achieving the property is generally difficult because the necessary and sufficient 
conditions are restrictive. Furthermore, there is a negative result that the non- 
AS-origin never becomes the UASAS-origin by the addition of any diffusion 
term m-, that is, as long as we employ the stabilization by noise with the use 
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of Wiener processes, we should allow the possibility of all sample paths traveling 
to points very far away from the origin. 

On the other hand, we claim that the stabilization by noise has still possibil¬ 
ity to provide the non-AS origins with AS, provided that the resulting systems 
are neither ordinary nor stochastic differential equations. To achieve the result, 
we consider the addition of deterministic noises having unbounded variations. 
This plan needs to employ a system representation by using rough path analysis 
[3 m [13 [S] because the theory enable us to consider the dynamics including 
signals having unbounded variations such as Wiener processes and a particu¬ 
lar kind of deterministic processes. The key point of the analysis is to classify 
external inputs by calculating the orders of the variations that is finite. Based 
on the information of the orders, rough paths and their dynamics—generally 
said to be rough differential equations, and rough systems in this paper—are 
obtained. 

In this paper, we provide a concrete system representation of rough systems 
in show that asymptotic stability in roughness (ASiR )—this is the property 
for the origins of rough systems—is the same property as UASAS, and demon¬ 
strate that ASiR is available for stabilization strategies by using deterministic 
signals having unbounded variations. Especially, we confirm that there is a case 
that such deterministic noises are capable of stabilizing the origin of dynamical 
systems in the sense of ASiR while stochastic noises such as Wiener processes 
do not achieve the origin being UASAS. 

The rest of this paper is organized as follows. In Section |2l we briefly state 
the motivation of our research by showing the strategy of stabilization by noise 
using stochastic and deterministic control inputs. In Section |3l we summarize 
the previous results of rough path analysis; especially we show the definitions 
of geometric rough paths, integration along rough paths, and rough systems. 
Then, Section |4] shows main results of this paper; in Subsection 14.11 we derive a 
concrete system representation for our purpose by using canonical rough paths; 
in Subsection 14.21 we provide ASiR property; in Subsection l4.31 we confirm that 
the ASiR is the same property as UASAS; furthermore, in Subsection 14.41 we il¬ 
lustrate the answer to our motivational example and discuss a further advantage 
of the usage of deterministic noises against stochastic ones. Section |3 concludes 
this paper with some important remarks for control problems. 

Notations. Let M" denote an n-dimensional Euclidean space; in particular, 
R denote If 7 : [0,oo) —>■ [0, 00 ) satisfies 7 G AToo, then 7 (r) is monotone 
increasing, 7 ( 0 ) = 0, and 7 (r) —)■ 00 as r —)■ 00 . Further, wt G is a d- 
dimensional, independent, and standard Wiener process; i.e., the values of the 
variances are all t, and all the covariances are constantly zero. The differential 
forms of ltd and Stratonovich integrals of a : R” R^^ in Wt G R'^ are denoted 
by a{x)dw and a{x) o dw, respectively. For w : R" —>• R, gi, 32 : R" ^ R” and 
a; G R", we use the notation of {Lg^v){x) = {dv/dx){x)gi{x) and {ljg.^hg.^^v){x) = 
{dljg^v/dx){x)g2{x). 

For the usage of rough path analysis, let J = [0, T] with T > 0; a subdivision 
of J is an increasing sequence of real numbers D = {to,ti,... ,tN) with N = 
1 , 2 ,... such that 0 < to < U < ••• < < T; and \Dt\ = max{to,ti — 
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to,... ,t]\j — iw-i ,T — Tjv} is a maximum value of lengths of time subintervals 
for some fixed -D; further, 

At := {(s, t) e [0, T]^\0<s<t<T} (1) 

be a set of a pair of the initial and the terminal values for time variables; and 
Aoo is At with T oo. Let also 

T2(]R") = K0R" © (K" OR"); (2) 

that is, an element of T^(]R") consists of a scalar, a vector, and a matrix. For 
any A, B £ T^(]R") with n = 1,2,... and A = 0,1, 2, a product of A and B is 
defined by A ® S = C G T^(]R"), where 

A 

:= ^ A'® (3) 

1=0 

For any spaces M and TV, the coordinate projections are defined by ttm ■ M © 
N ^ M and ttm : M (B N ^ N. For any function yr : [0, oo) Y for any set 
F, j/s.t : At —>■ F is defined by ys,t = Ut — Us- The i-th element of a vector 
a G is represented by a\i\, and the t-th row and the j-th column element of 
a matrix A G R'^ © R^^ is denoted by A[i, j]. For c G R, 5° means b to the power 
of c if & G R, and denotes the c-th element of i? if i? G r^(R"’) (see also 
Definition [5]). 


2 Motivational Example 


To begin with, let us consider a stabilization problem for the origin of the 
following simple system: 

it = fixt) + gi{xt)ut[l] + g2ixT)ut[2], (4) 


where 


fix) 


-7 

0 


^ X, gi{x) 


0 

1 


0 

0 


X, g2{x) 


0 1 
-4 0 


(5) 


Ut G R^ is a control input vector, xt G R^ is a state vector, and t G [0, oo). 
The origin of the system can be stabilized by just the addition of noises. For 
example, if we employ u as a Gaussian white noise, then (|3]) becomel^ 


dxt 


f{xt)dt + gi{xt) o dwt[l] + g 2 ixt) o dwt[2] 


9 

0 


dt © 


0 

0 


xtdwt[l] + 


0 

-4 


1 

0 


xtdwt [ 2 ]. 


( 6 ) 


^More precisely, we design u = , where G is defined in Definition!?] then we 

obtain lO by the limiting operation as with Theorem [5] 
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Analyzing the system using the basic results of stochastic Lyapunov stability 
theory, all the sample paths of this system converges on the origin with proba¬ 
bility on^. Fig. [1] shows an example of the sample paths. 

Emphasizing again that the unstable origin becomes the stable origin by just 
adding the Gaussian white noises—roughly speaking, the noises are “deriva¬ 
tives” of Wiener processes. There are two candidates for the cause of the stabi¬ 
lization result: the randomness of the Wiener process and the unboundedness of 
the variations of the processes. In this paper, we claim that the main cause is the 
unboundedness, generally not the randomness because there are deterministic 
control inputs u = u^{r]), where 


Ut [!](??) = b 


cos(77^f) — 1 


[ 2 ](??) = b2 


sm{r]^t) 


(7) 


with T] = 1 , 2 ,..., deriving a stabilization result similar to the above example; 
see Fig. [5] as the circumstance evidence. As r] gets larger, the trajectories seem 
to close with a particular smooth trajectory that converges on the origin; to 
investigate the true identity of the “pseudo-smooth trajectory”, we need consider 
the situation of ry —>■ oo. However, if ry —>■ oo, the derivatives of ^^[l] and 
are impossible to be defined because they become signals having unbounded 
variations, as with Wiener processes. 

Thus, to deal with “deterministic unbounded-variation noises” such as lim^^co u{r]) 
as with “stochastic unbounded-variation noises” such as Wiener processes, a 
comprehensive representation method of deterministic and stochastic systems 
is needed. Fortunately, the recent works on rough path analysis [nmiiiT] en¬ 
able us to deal with such unbounded-variation noises in a unified way. Then, 
in m. Lyapunov-like stability theory for rough path analysis, namely ASiR, 
has been proposed. In this paper, we develop the theory by providing the 
proofs for the related rough systems and making further discussions so that the 
analysis of stabilization by noise is more generalized by allowing deterministic 
unbounded-variation noises. Furthermore, we will achieve that deterministic 
unbounded-variation noises provide stability property tighter than stochastic 
unbounded-variation noises. 


3 Rough Path Analysis 

To achieve our system representations including both deterministic and stochas¬ 
tic unbounded-variation noises, we should consider state variables and their dy¬ 
namics allowing such noises. Because we claim that rough paths are suited for 
our state variables, we briefly summarize the rough path analysis based on [l6] 
in this section. 

^The origin is globally asymptotically stable in probability because Vi(x) = x^x is a global 
stochastic Lyapunov function; briefly speaking, {CV-i){x )—as the same form as II33I I below—is 
equal to —(l/2)x^x; it is negative definite for all x S R^; see the references on stochastic 
stability, for example n Emu US]. 
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Figure 1: A sample path of Q with u = 
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(b) r] = 10 



(c) 77 = 100 

Figure 2: Paths of (HJ with u = [rf) , bi = 2> and 62 = 4. 
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3.1 Rough Path 

Here we start with the definition of p-variations that are necessary to achieve 
the notion of rough paths: 


Definition 1 (p-variation (16) i Let p > 1, a continuous path x : Ay —K", 
and a sequence of time sub-intervals D be considered. Then, 


\p,J ~ 


N-l 

sup V \xt^,tk+i\^ 


1/p 


( 8 ) 


is said to be the p-variation of x (on the interval J). Furthermore, */||a;||p_j < 
00 , then x is said to have finite p-variation (on the interval J). □ 


To simplify the discussion, we consider the situation of p S [2, 3) in what 
follow^; this assumption also reduces the complexity of definitions hereafter. 
The next definition describes the predecessor of rough paths: 

Definition 2 (multiplicative functional (16) 1 Let X : At r^(]R") he a 
continuous map satisfying 

Xs,t = {Xl,Xl„Xl,) (9) 

for each {s,t) G At, where X°^ G K, X^^ G K”, and Xf ^ G K" ® K". The 
function X is called a multiplicative functional (of degree 2 in K"J if ^ = 1 
for all (s, t) G At and 

Xs,T ®Xr,t = Xs,t, Vr G [s,t]. (10) 


□ 


Now we are ready to define rough paths as follows: 

Definition 3 (rough path (16) 1 Let X, X{1),..., X{ri) : At — t T^(R”) be 
continuous maps. Then, 

1. Xs^t is said to be a rough path (of degree p in if it is a multiplicative 
functional with X^ having finite p-variation; 

2. A(l)s_t,..., A(77)s_t are said to be smooth rough paths if they are multi¬ 
plicative functionals with A(l)^,..., X{r]f having finite 1-variations, re¬ 
spectively; and 

®Of course, we can consider rough path analysis for any large p > 1; however, considering 
p € [2, 3) is sufficient to discuss the relationship between deterministic noises ([7jl and Wiener 
processes constructed by Theorem [5] 
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3. Xs^t is said to be a geometric rough path if it satisfies dp{X{r])s^t, Xg^t) —^ 0 
as r] ^ oo, wher^ 



V 


( 11 ) 


Furthermore, an element Xf ^ for j 


0 ,1, 2 is said to he the j-th level path of 


□ 


Note that, we mainly consider the first level paths of geometric rough paths 
as state variables; rough paths and smooth rough paths are introduced for the 
definition of geometric rough paths; and the second level paths are needed for 
describing the dynamics of the first level paths; see also the discussions in the 
rest of this section. 

Remark 1 The second level paths of rough paths generally have flexibility of 
the definitions; for example, if X^ ^ = w{t) — w{s), Xf ^ can be defined by 
using whether ltd integral X^ ^ = fgW(T) 0 dwfr) or Stratonovich integral 
Xf t = fg u’(t) 0 odw(T). However, if Xf j. are created by ltd integrals, X 
is not a geometric rough path because there is no smooth rough path satisfying 
dp(X(ri)s,t,Xs,t) 0 as rj ^ oo, see Sec. 3.3 in In contrast, if X'^ is 

constructed by Stratonovich integrals, X is a geometric rough path because the 
approximation theorems |21 17^ typified in Wong and Zakai ]25j . imply that 
the existence of X{r]) such that dp{X{r])s^t, Xg^t) —t 0 as rj ^ oo; see Sec. 3.4 in 
fW/ or l^. 

3.2 Integration Along Rough Paths 

Because we will consider (the first level paths of geometric) rough paths as 
state variables of dynamical systems, we should investigate the behavior of 
rough paths concerned with time. This aim needs differential equations for 
rough paths; however, the equations are incapable of being defined directly be¬ 
cause the first level paths of rough paths generally have unbounded 1 -variations; 
that is, they are undifferentiable almost everywhere, such as Wiener processes. 
Therefore, we should consider to derive integrals along rough paths and their 
integral equations as with stochastic systems. According to [16], we show the 
definition of integrals along rough paths. For s = to < ■ ■ ■ < tk = t, {s,t) £ At 
and a smooth function h = {hi, /i 2 ,..., hq) with hi, / 12 ,..., hg : K.'J —>■ R'J, let us 
consider 

4To be precise, we have to define p-variation topology in Definition 3.12 of m after the 
definition of control functions, which is a completely different notion from control inputs in 
this paper, in Definition 1.9 of m- While they are necessary to provide geometric rough 
paths and prove Theorem^below, we abbreviate them to simplify the structure of this paper. 


7 



Ih{X,D)[y. 

N 

k^l 

(12) 

with Z = 1,2 and 




■.=hYxl_Yxl_^,^ + {VH){xl_^^,j, 
■.=h{xl_Y®h{xl_Y-xl_^,^, 

(13) 

(14) 

where 




Fill ■ 

JW:=tr 

(15) 


for i = 1, 2,..., q. Then, we define the following: 

Definition 4 Let Xg^t ■ T^(]R") be a geometric rough path and h{X^) be 

smooth enough. Then, 


[ h{Xr)dX\:= lim Ih{X,D)it 

Js IDtKO 

is said to be an integral of h along X^ for I = 1,2. Further, 
h{Xr)dXr-.= (^1,J^ h{Xr)dX^,J^ h{Xr)dxi^ 
is said to be an integral of h along X. 


(16) 


(17) 

□ 


Remark 2 The key point of the definition along rough paths is the usage of 
(fT^ .- Ih{XY is considered including the second order Taylor polynomial of 
h{X^) if we consider X^ £ K, Xl ^ = Xf — Xj and X^ = {Xf — XlY/2. 
Roughly speaking, the integration along rough paths is an extended notion of 
integrations by considering higher-order Taylor expansions. ♦ 


3.3 Rough Systems 

Through the definitions of rough paths and integration along rough paths, we 
eventually achieve representations of dynamics for rough paths. Let us consider 
a nonlinear systenH 

N 

XB,t = ( 18 ) 

fc=i 

where {s,t) G At, x := (a;[l],... ,x[n])^ : At Mx = K” is a state vector, 
u := (m[ 0 ],u[ 1 ], ... ,u[m])'^ : At —>■ Mu = is an input vector, ^,-[ 0 ] = r G 

^If u has finite l-variation, IjlSIl is equal to i = g{x)u. 
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[s,t], and g = {go,gi, ■■■,gd) with go,gi, ■■■,gd- Mx Mx are assumed to be 
smooth and locally Lipschitz. Furthermore, we only consider the combination 
of g and u providing a unique global solution to mB 

Let a geometric rough path U : At —>• T‘^{Mu), a map X : At — >■ T‘^{Mx), 
and a map Z : At —>• T^{Mu)®T^{Mx) satisfying = U, ttmx{Z) = X; 

i.e., 


Kt = 


ult 

Kt 


Zlt = 


ult 

vXU 


Kt 


(19) 


be considered for s < t < t e [0,T] and : At —>■ M" 0 K", : At 

R™ (g) R”, and : At -)■ R" 0 R™. 


Definition 5 (rough system and rough solution) An equation 


Zs,t — 


9iK) 


dZr 


( 20 ) 


and Zsd o,re said to be a whole rough system and a whole rough solution, of 
dni), respectively. Furthermore, 


X 


1 _ 


ttmx 




( 21 ) 


extracted from (1201) and X^ are said to be a rough system and a rough solution, 
of (US, respectively. □ 


The following result is a simplified version of Lyons universal limit theorem: 

Theorem 1 ( |16] 1 Let u(rf) for g = 1,2,... have finite 1-variations and U(rf) 
be a smooth rough path satisfying U^{g)s,t = Usdiv)- Let also Us,t be a geometric 
rough path satisfying dp(Usdiv),U) 0 as g ^ oo and the initial value Zg G 
T2(M[/)©tLpJ)(Mx) be fixed. If dTS]) has a unique global solution with u = u{g) 
for any g, the whole rough solution Zg^ of m uniquely exists and Zg^t is a 
geometric rough path. ♦ 

This result immediately implies that X € T^(Mx) is also a geometric p- 
rough path of x. In this way, we obtain m as the dynamics of our state 
variable note that, Theorem [1] implies that we should define the whole 
rough system (uni) because the dynamics of X^ is directly affected by see 
also Remark 2] below. This fact will be clearly articulated by deriving more 
specific representations in the next section. 

®In 1161 . g is assumed to be globally bounded for the existence and the uniqueness of global 
solutions. However, we consider g satisfying local bounded conditions; the lessening of the 
conditions on g is necessary to consider our problem formulation, for example, (O in Section [2] 
The detailed discussion on more concrete assumptions on g and u will be described in the 
next section; see B1 and B2 below. 
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Remark 3 Note that, in HEW, a whole rough system is said to be a rough 
differential eguation with the abbreviation of integral signs “f ” and X is said to 
be a rough solution. ♦ 

Remark 4 Despite focusing attention on the dynamics of , we have to define 
the whole rough system (Eni). The reason is that J XdU is incapable of being 
defined if X and U have different ranges, i.e., T‘^{Mu) 7 ^ T^{Mx) while f ZdZ 
can be defined UMW- Furthermore, The ltd map of X^, i.e., 1 —7 X^ (note 

that J7t[0] = t), is generally discontinuous while one of Z, i.e., U 1 -^ Z, is 
continuous —the latter map is said to be Ito-Lyons map M- In theory, 

this point is the true identity of “the bridge” between deterministic and stochastic 
differential equations. ♦ 


4 Stabilization by Unbounded Noises 

In this section, we investigate the effects of unbounded-variation noises—such as 
the limitation of u{ri) as ry —>■ 00 in Section[5]and Wiener processes—on stability 
properties for dynamical systems. 


4.1 System Representation 


To begin with, we formulate rough systems more specific than the previous 
section. 

In what follows, Let us consider the followings: 

AO ft( 77 ) has finite 1 -variations for all ry = 1 , 2 ,...; 

A1 U{r]) : Aoo —>■ T^{Mu) is a smooth rough path satisfying 


Ultiv) =Us,tiT), Ultip) 


Ms.r(?y) ® du{ri)r 


( 22 ) 


for all ?y = 1 , 2 ,...; and 

A2 Us,t '■ Aoo —>• isageometric rough pathsuch that dp([/(?y)s^t,[/s_t) —>■ 

0 as ry —>• 00. 

Then, the rough system is derived as follows: 

Theorem 2 Let A0-A2 be assumed and u = u{r]). Then, the rough system of 
(HHi) is given by 


/ t m m m ^ 

Y,9j{Z()dU([j]+ / Y.J2^(^r)9k{X()dU^[k,j] 

j=0 


as long as it has a unique global solution. 


(23) 

♦ 
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The proof is shown in Section|B] note that, although the assumptions A0-A2 
have already appeared in Section 3 of m, the concrete system representation 
shown in this theorem was not provided. The rough paths satisfying AO are 
said to be canonical rough paths [S], and they are used for dealing with Wiener 
processes in rough path analysis; see also Subsection 14.31 In this paper, we also 
consider AO for constructing ASiR properties. 

Remark 5 Emphasizing again that, not only the rough solution is uniquely 
derived but also the rough system (1231) includes the second level variations: there 
exist the second term of the right-hand side of (ED; of course, if u has finite 
1-variation, the terms vanish and (1231) is equal to the original system dUD as 
usual sense of Stieltjes integrals. Furthermore, if u has finite p-variation with 
p > 3, we add terms of the third level variations to the right-hand side of (I23p 
based on (Uni- ♦ 

4.2 Lyapunov Stability for Rough Systems 

In this section, we consider stability properties for the origin of the rough system 
(1^ : that is, we consider po(0) = 0 and a control input u satisfying A0-A2 and 
the following: 

A3 Us,t makes the origin of (1231) a unique equilibrium. 

Here we define the notions of stability for the origins of rough systems. 

Definition 6 (stability in roughness) Let Us,t be fixed such that AO-A3 hold 
and Zg be fixed for an initial time s > 0. The origin of (US with u = 
hm^^oou(p), or (l23l) . is said to he 

El(6]-1 (uniformly) stable in roughness if for each e > 0, there is S = d{e) > 0 
such that 

|x1|<^=^|a1|<£ (24) 

holds for all t > 0, uniformly in s > 0 ; 

El [^2 (uniformly) locally asymptotically stable in roughness (locally ASiR) if it 
is stable in roughness and 6 can be chosen such that 

\Xl\<6^ \im Xl = 0 (25) 

t—yoo 

holds uniformly in s > 0; 

El(6]-3 (uniformly) globally asymptotically stable in roughness (globally ASiR) if 
it is stable in roughness and moreover 

lim X] = 0 (26) 

t—^OC 

holds for all X) G R", uniformly in s > 0. □ 


11 


This definition is the same as Lyapunov stability, local AS, and global AS 
for the origins of ordinary differential equations 0110 ]; hence, a question arises 
whether the related Lyapunov stability theory for Definition[0]is also the same or 
not. To investigate it, we consider the dynamics of a sclar function v : K" —>• R. 
as a candidate of “Lyapunov functions”; let us consider the following enhanced 
system: 


N 


Xs,t — 'y ^ 

k=l 

where x = (a;[ 0 ],a;[l],... ,x[n])'^ and 
g[x) = 


{hg„v){x) ... {l^g^v){x) 

9{x) 


(27) 


(28) 


Thus, we obtain the following: 

Lemma 1 Let A0-A2 he assumed. Then, we obtain 

nt TYi pi m m 


pi, I I 0 I I 0 I I L 

Kt= / ^(L,,l/i)(A,i)dt/i[j]+ / ^^(L,,L,^l/i)(A;)dC/2[fc,j], 

S _1 S ^ _n 1 _r> 

(29) 




j—0 k—{) 


where := A^[0] is said to be the fist-level path of a rough path V, as long as 
exists as a unique global rough solution of (EZj). ♦ 

The lemma is proved as with Theorem 0 because the form of is the 
same as m, the rough system is o with replacing Z by Z = (U,X) and g 
by g, respectively. Thus, we obtain (1^ . 

Now we assume the following with D as an open subset of R" including the 
origin a: = 0 : 

Cl V : D ^ M. is smooth for all x G D and there exist continuous, positive 
definite, and proper functions ITi, W 2 : D —>■ R such that 


Wi{xt) < v{xt) < W 2 {xt) 


(30) 


for alH > 0 and x G D. 

Then, we obtain the following result: 

Theorem 3 Let A0-A3 and Cl he assumed; let also dUD and (1291) be considered 
with Xs being fixed for an initial time s > 0. Moreover, assume that there exists 
a function DV^ : R" —)• R such that = f*('DV^)(Xf)dT holds in D. Then, 
the following results are obtained: 

C|3|-l if {'DV^){X^) < 0 for all X^ G D and t > 0, the origin of (1181) is stable 
in roughness; 
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C|^2 if there exists a continuous positive definite function W3 : Z? —> R such 
that {'DV^){Xf) = —WalXf) holds for all X^ G D and t > 0, the origin 
of ® is locally ASiR; 

Cl^3 if Cl and (J^2 are all satisfied for D = R", the origin of (I18p is globally 
ASiR. 4 

The proof is the same as one for the original Lyapunov stability, for exam¬ 
ple, Theorems 4.1 and 4.2 in m because = f*{'DV^){Xf)dT is the same 
condition as Lyapunov functions for ordinary differential equations. Thus, we 
obtain stability properties for rough systems; we conclude that the function v 
of Cl, or the first level path of the rough path V, is said to be a rough 
Lyapunov function (RLF), a (strict) local RLF, or a (strict) global RLF, if it 
satisfies CHI, 2, or 3, respectively. 

Remark 6 Assumption AS is natural because we consider the stability of a 
unigue equilibrium; however, this assumption does not ensure that the origin 
of dm with u{r]) is an equilibrium when rj is finite. This point makes strange 
phenomena on the stability analysis of rough systems; see Subsection \4-.4\ 4 

4.3 Relationship with stochastic stability 

The previous subsection shows that stability in roughness is the same notion 
as stability for ordinary differential equations, provided that the system are not 
described by ordinary differential equations but rough systems. On the other 
hand, rough systems include stochastic systems; hence, stability in roughness 
should also be the same as some notion of stochastic stability. In this subsection, 
we investigate what is such stochastic stability. 

In this subsection, we assume that u{ri) and (fTSll satisfy the following: 

B1 The signal u(rf) is piecewise linear for all time intervals. Furthermore, if 
b] = with any constant ai,..., a™ G K for any j = 1,2 ... ,m 
and for any time interval, (I18|) has a unique global solution. 

Lemma 2 Let (1151) and (1^ be considered. Let u = [rf), where 

is defined in Definition^ If B1 holds and an {n1)-dimensional stochastic 

system 


Xs.t = 


pt rn 

/ 9o{x)dT + '^ gm{x)odWr[j] 

J S J S 


has an unique global solution, then we obtain 

pt "1 pt 

Vl,= / {CV^){Xl)dT+ Y, / [Lg^V^){Xl)dUl[j\, 

J S X_1 s 


t=l • 


(31) 


(32) 
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where 


{CV^){Xl) = iLgX)iXl) + (33) 

i=i 

♦ 

The proof is based on the discussion of Section 3 in [TC] ; see Section [C] below. 

This lemma immediately means that is the same as the infinitesimal 
operation of for Stratonovich-type stochastic differential equations driven by 
standard Wiener processes m- Thus, we obtain the following: 

Proposition 1 Let dUD and (l29)) be considered. Let also u{r]) be as with 

Lemma\^ If B1 holds and 

{Lg^V^){Xl)=0 (34) 

holds for all j = 1,2,... ,m and for all X^ G K" and t > 0, then we obtain (1331) 
andCV^=VV\ ♦ 

Because uniform almost sure asymptotic stability (UASAS) [3] requires con¬ 
ditions as with (IMl) . this corollary implies that UASAS is the same property as 
ASiR. It is not a surprising result because the UASAS is “almost the same” as 
AS [TO] . 

At the same time, in |19j . it is proved that the addition of any Wiener pro¬ 
cess never makes the non-AS-origin become UASAS. In contrast, we claim that 
deterministic unbounded-variation noises, such as are capable of stabilizing 
the origin in the sense of ASiR, despite the stability being the same as UASAS. 
The next subsection shows positive effects of deterministic unbounded-variation 
noises by comparing with stochastic unbounded-variation noises. 

Remark 7 If B1 holds, there exists a unique global solution for under 
u = w^^{r]) with finite rj, where w^^{r]) is “piecewise linear approximate 
Wiener processes”; see Section HI below. Further, the global solution of dSD 
is needed for the existence and uniqueness of ([32]). The existence of the both 
unique global solutions enable us to measure dp{X{rf),X) as rj ^ oo; thus, we 
are capable of using Theorem [I] Note that, B1 allows us to consider bilinear 
systems such as o because, if all the control inputs m[1], ... ,u[m] are nothing 
but step inputs, then o becomes a linear system; that is, the system never 
has any finite escape time for any initial state. Of course, if we employ more 
complicated construction procedure for Wiener processes, B1 is insufficient to 
ensure the existence and uniqueness of the solutions. ♦ 

Remark 8 As stated in Section there are other schemes for constructing 
Wiener processes not satisfying Bl. Nevertheless, the results in m imply that 
no Wiener process make the non-AS-origin become UASAS; that is, as long as 
considering the negative result, it is enough to consider Wiener processes created 
in Section\^ ♦ 
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4.4 Stabilization by Deterministic Unbonnded Noises 

Here we describe stabilization results by the addition of control inputs having 
unbounded variations with reconsidering the discussion on Section [21 In this 
subsection, we assume that u{ri) and (IT^ satisfy the following: 

B2 The signal u{f]) satisfies {v)[j] —>■ 0 as 77 —>■ 0 for any j = 1,2 ... ,m 

and for all time intervals. Furthermore, if = 0 for any time inter¬ 

val, (ITSl) has a unique global solution. 

The following theorem clarifies “the additional terms” generated by the de¬ 
terministic unbounded-variation noises 0: 

Theorem 4 Let us consider (HID with m = 2 and ©. Let also assume that B2 
holds and the integral form of an n-dimensional ordinary differential equation 



has a global solution. Then, the rough system is derived by substituting x = 
into (I35]) . ♦ 

The proof is shown in Section [P] below. 

Remark 9 If B2 holds, there exist unique global solutions for (HI under u = 
u{r]) with finite rj. For example, let us consider n = 1, m = 1 and g{x) = (0, x^); 
then we have x = x'^ulrf); that is, the solution is Xt = Xs/{1 — XsUs,t{v))- This 
implies that there exists a solution for all timet € [0,oo) if\us,t{'n)\ < l^;®! holds; 
this inequality always holds with sufficiently-large rj because B2 implies that the 
amplitudes of uljf) comes close to zero as p increases. This implies that, there 
exists r]' = t] — a with some positive integer a such that u{rj' -\- a) ensures global 
solutions for all rj' = 1,2,.... On the other hand, the global solution to (I55|) 
is needed for the existence and uniqueness of the related rough system. Note 
that, the existence of global solution to dSSl) is ensured if g is locally Lipschitz 
and there exists a global Lyapunov function; if finite escape time arises for some 
initial state, it contradicts the existence of a global Lyapunov function. In this 
way, we can measure dp{X{r]), X) as rj ^ 00 to employ Theorem[l[ ♦ 

4.5 Discussion With Case Study 

Here we reconsider the motivational example of deterministic noises in Section |21 
Using Theorem m we obtain 

yi _ ^ —14-I-6162 0 

^’*"27, [ 0 2-6162 

as the rough system of ([3]) with u = u^{rj). Therefore, if 6162 S (2,14), the 
origin is obviously globally ASiR because V^[X^) = (1/2)(X^)^A^ immediately 


Xldr 


(36) 
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0 0.2 0.4 0.6 0.8 1 1.2 

Time 


Figure 3: A path of with u = u^(r/), rj = 100, = 3 and 62 = 4. 


satisfies the condition C[3]-3. Fig. [3] describes the trajectory of this for the 
confirmation. 

The answer to the motivational example shows that deterministic unbounded- 
variation noises are capable of making the non-AS-origin become the ASiR- 
origin. This implies that, the example is the case that stochastic and deter¬ 
ministic unbounded-variation noises both available for stabilization by noise; 
however, the results of m imply that stochastic unbounded noises do not help 
the non-AS-origin become the UASAS-origin. Thus, deterministic unbounded 
noises have possibility to achieve the stabilization results stronger than stochas¬ 
tic unbounded noises. In this subsection, we consider the reason by using a 
simpler example. 

Considering n = 1, 61 = 62 = 1, and 

go = 0, gi{x) = 1, g 2 {x) = -x^, (37) 


the rough system (1551) is calculated as 


Ai, = - [ Xldr 

J S 


(38) 


while f* dU^ = 0. Obviously, the origin is globally ASiP. 

Of course, the trajectories vibrate because Uiig) and are the deter¬ 

ministic noises. However, the origin of the related rough system is globally 
ASiR. More clearly, if g is finite, the system is represented by 


^s,t — 


= / u^[Mig)dT- 


f 


xiu^ 


[2]{g)d7 


(39) 


that is, as long as g is finite, the origin is never an equilibrium; of course, it 
is not globally AS. Nevertheless, the rough system (1551) implies that the origin 
is globally ASiR. The trick is made by the calculation results ^ 0. Thus, 
the representation of rough systems enlarges dynamical systems that we can 
consider. 
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The above stabilization result is due to the form of Theorem 2] for compar¬ 
ison, if M = is chosen for with dSil), the resulting stochastic system 

is 


f odwr[l] — f a;^ o dwi-[2], 

' s J S 

j-t j-t j-t 

(40) 

/ dwT-[l] — / x^dwrl^] + x^dr 

' s J s J s 

(41) 


if a solution existfl Because the origin is not an equilibrium due to the form 
of the first term of the right-hand side of (HD), it is not stable in probability. 
That is, the addition of Wiener processes makes the resulting stochastic systems 
have diffusion terms, namely the first and second terms of the right-hand side 
of (HTl) . On the other hand, the addition of deterministic noises ([T]) does not 
cause such terms in rough systems, see again (l35l) . This is the true identity of 
the advantage of deterministic unbounded-variation noises lim,j_>oo {v) over 
stochastic unbounded-variation noises w. 


5 Concluding Remarks 

In this paper, we provided the availability of “deterministic unbounded-variation 
noises” by describing rough systems and their Lyapunov stability properties. In 
the procedure, we also clarified that the noises are capable of ensuring asymp¬ 
totic stability in roughness, which is generally a stronger property than asymp¬ 
totic stability in probability. Furthermore, the noises sometimes stabilize the 
origin of rough systems while the original system with bounded noises—they 
are considered as approximations of deterministic unbounded-variation noises— 
does not have the origins in equilibrium. 

Finally, we should notice that the analyzed issue is similar to the problem 
formulation of stabilization by the approximation algorithms OHS]. The al¬ 
gorithm has not clarified the relationship between deterministic and stochastic 
noises yet. The comparison should be important to confirm that our strategy is 
effective for the application of stabilization by noise to nonholonomic systems 
[20] . This point will be solve by considering rough paths of higher-order degrees 
with some improvement of the formulations. 
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Wiener Processes Created by Wong and Za¬ 
kai |25| 


In this section, we consider a creation method of Wiener processes introduced 
by Wong and Zakai pS] . 

Definition 7 ( (25) ) Let N = rj. Let also bring on sequences of discrete-time 
Wiener processes {wP^[j]{T]), ..., wP^[j]{r])} for j = 1,2,... ,m and 4 G (s,t) G 
At. For every r G [tk,tk+i), 

(42) 

— tk 

Then, w^^[j]{ri) for j = 1,2,..., m is said to be a (approximate) Wiener pro¬ 
cess. □ 


Theorem 5 ([25]) Let m = 1,2,.... If u = w^^{'q), (IT8)) becomes 

rt 2 


B,t= / go{Xr)dT / gj{Xr) odwrij] 

J S 

= / fixT)dT + '^ / gj{Xr)dWr[j] 

S X_1 ^ s 


(43) 
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as T] ^ oo with probability one, where 


f{Xl) := g,{Xl) + ^ --^^^xl)g,{Xl). (44) 

i=i 

♦ 

Remark 10 The above-mentioned Wiener processes are piecewise linear in¬ 
terpolations of discrete-time Wiener processes; that is, satisfies B1 while 

rj = N is finite. Note that, there are some creation methods of Wiener processes, 
see IB [MW- If an approximate Wiener process is constructed by nonlinear 
interpolation functions, then the derivative does not consist of step inputs; in 
this case, B1 does not hold. ♦ 


B Proof of Theorem [2] 

Let Z : At T^{Mu) © T^{Mx) be the whole rough solution to (fT51) with 
•kmu{Z) = Uig) and for fol g = 1,2,.... In this section, let 

s' = tk-i and t' = tk for any fc = 1, 2,..., TV for simplisity. 

B.l Deriving a geometric rough path Z. 

First, we calculate the elements of the rough path Zs'^t' ■ The first level path is 
immediately obtained as Zl, = lim^^oo Zl, . The second level path is derived 
as follows. For i,l = 1,2,... ,n and j = 1, 2,..., m. 


Y^V^[J,^]= hm [\uig))l,^Aj]d{X{g))l[z] 

m 

= lim ^5fc[t]((l(r7))i,) / {U{g))], [j]d{U{g))l[k] 

m 

= lim '^gk\i]{{X{g))l,){U{g))l,^t'[Tk] 

■n—>-oo ^ ^ ’ 

m 

= Y^gk[i\{Xl)Ul^AjM- 


(45) 

(46) 

(47) 

(48) 


Note that, the above integrals in (l45ll - (l46l) are all Stieltjes integrals because 
U^{g) and X^{g) are nothing but the smooth functions u{g) and xlr;) that 
have finite 1-variations, respectively, where u{g) is assumed as AO, U^{g) is 
defined in Al, and x{g) is a solution to the original system (fTSl) with Us,t = 
u{g)s,t. Transforming (l47l) into (l4^ is achieved by A2 and the assumption of 
the existence of a unique global solution x{g) for any g > 0. 
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Similarly, we can calculate 

m 

[j, i]='^9k w {X},)U^,^t,[k,j], 

k^O 

m 

X!>A^,l]= 9jmxh)9miXh)UlAj,k]. 

j,k=0 

Thus, all the elements of Z^, are calculated. 


(49) 

(50) 


B.2 Deriving a rough system. 

Using the form of (fT^ with 

o). 

g = n + m + 1, and 2 = (m[0], ..., u[to], x[l],..., x[n])'^, we obtain 

m 

{h{Zl,)Zl,^,,M=Y,gMZs')Ul,^A^ 

3=0 


and 


m m+n 


(vi?)(x4_^,jw=EE 


j—O k—0 
m n 


dZ^[k] 


ixl,)zl,Ak]uyAj] 




dX} 

j=0 z=i ^ 
m m n 


= E E E ^Ml')9k[l\{Xl,)Ul^,[k,j\ 

j—0 k—0 1—1 ^ ^ 


(51) 


(52) 


(53) 


for i = m + 1,..., m + n. Thus, we obtain Ih{Z)\, j,. These calculation results 
and the definition of the integral along rough paths in Subsection 13.21 implies 
that (1^ holds. 


C Proof of Lemma 

If AO and A1 hold, all the elements of the second level path of U{rj) are trans¬ 
formed into 

Ult[hk]{T]) = ]^UsAj]{ri)us,t[k]{9) + As,t[hk], (54) 

where _), /c = 0,1 ,..., m and 

y UsAi\ijl)duT[k]{'n) -1^ j UsAk\{v)dur[j]{'n)- (55) 
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On the other hand, if we design Ur(v) = ,..., , we 

obtain 


lim wY^ijl) = , 

r]—¥co 

lim [ [j] {ri)dw^^ [j](T])=t-s, j = 1 , 2 ,..., m 

V^oo 

lim f [j]{r])dw^^ [k]{r]) =0, k^j, j = k = 0, 

by using Ito’s stochastic analysis. Therefore, we obtain 
= (t - s,Ws,t[l], ■ ■ ■ ,Ws,t['m]V ^ 

UltijJ] = - s), j = 1, 2 ,..., m 

Ult[j,k] = 0 , k^j, j = k = 0. 


(56) 

(57) 

(58) 

(59) 

(60) 
(61) 


with As^t[j, fc] = 0 for all j,k = 1,2,... ,m by using (1551) . Then, substituting 
these results into (15^)) . we obtain dSSl). 


D Proof of Theorem [4] 


To begin with, we consider u = (u[0 ],m[ 1],zi[2])^, where Ut[0] = t and u[j] = 
u^{il)[j] for j = 1,2. As rj ^ oo, p-variation of u is infinite if p G (0,2) and 
finite if p > 2; that is, AO holds. Hence, we can consider Al, and then we obtain 
a geometric rough path Ug^t satisfying A2. 

Next, we calculate elements of U: 


UlW = UltM = UlAO, 0] = c/2 JO, 1] = c/2 JO, 2] 
= C/2J1,1] = C/2J2,2]=0, 

C/2J1,2]= lim (17(r;))2 [1,2] 

rj—^oo ’ 


(62) 


= lim / {u^[l]iri)-Ug [l]{ri)}u.,[2]{r])dT 


r]—^oc 


= lim / bib 2 {cos^{ri‘^T) — cos{r]‘^s) cos{'r]^T)}dT 


^^-^(t-s) = -UlA2,l\. 


Thus, we obtain 


(63) 


/ 

t — s 



0 

0 

o' 

Us,t = 1, 

0 


s) 

0 

0 

1 

V 

0 

Z 


0 

-1 

0 


( 64 ) 


Thus, we conclude that (1551) holds by using (IMl) and Theorem [5] 
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